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Airborne spectral measurements of surface-atmosphere anisotropy 

for arctic sea ice and tundra 

G. Thomas Arnold (SM&A Corp), Si-Chee Tsay (GSFC), Michael D. King (GSFC), 
Jason Y. Li (SM&A Corp), and Peter F. Soulen (UMBC/JCET) 

Popular Summary 

One of the factors affecting the Earth's climate is the reflection of sunlight 
from the Earth's surface back into space. The only practical way to measure 
reflected sunlight on a global scale is from satellite-based sensors. However such 
sensors typically measure the reflectance at only a single angle (usually from the 
surface directly below the satellite). This may or may not be representative of the 
true reflectance of the surface. To compute the true reflectance of a surface, the 
reflected sunlight from it in all directions must be known. The complete pattern 
of reflected sunlight from a surface is known as the bidirectional distribution 
function (BRDF). The more uniformly a surface reflects sunlight, the more 
representative the measured reflectance (such as by a satellite) is to the true 
reflectance. Thus for surfaces that do not reflect uniformly, knowing the BRDF 
for that surface is critical to knowing how to use the satellite sensor measurement 
to determine the true reflectance. In this paper we examine four common Arctic 
surfaces, snow-covered sea ice, melt-season sea ice, snow-covered tundra, and 
snow free tundra to determine their complete reflectance patterns and to 
determine how well a single reflectance measurement directly above the surface 
compares to the true reflectance of each surface. The reflectance measurements 
were obtained using an aircraft-based instrument called Cloud Absorption 
Radiometer during April 1992 and June 1995, and were done for several 
wavelengths of reflected sunlight. Polar plots of each measured wavelength 
showing the complete reflectance pattern of each surface are presented. The true 
reflectance for each surface is computed from this data and compared to the 
reflectance measurement directly above the surface. Reflectance above the snow 
surfaces most closely approximated the true reflectance, especially at the shortest 
wavelengths. Reflectance above the melt-season sea ice also fairly closely 
approximated the true reflectance, but slightly less so than the snow-covered 
surfaces. The snow-free tundra however was much different. Reflectance 
directly above the surface was as much as 50% less than the true reflectance. 



Abstract 


Angular distributions of spectral reflectance for four common arctic surfaces: 
snow-covered sea ice, melt-season sea ice, snow-covered tundra, and tundra 
shortly after snowmelt were measured using an aircraft based, high angular 
resolution (1-degree) multispectral radiometer. Results indicate bidirectional 
reflectance is higher for snow-covered sea ice than melt-season sea ice at all 
wavelengths between 0.47 and 2.3 pm, with the difference increasing with 
wavelength. Bidirectional reflectance of snow-covered tundra is higher than for 
snow-free tundra for measurements less than 1.64 pm, with the difference de- 
creasing with wavelength. Bidirectional reflectance patterns of all measured sur- 
faces show maximum reflectance in the forward scattering direction of the prin- 
cipal plane, with identifiable specular reflection for the melt-season sea ice and 
snow-free tundra cases. The snow-free tundra had the most significant back- 
scatter, and the melt-season sea ice the least. For sea ice, bidirectional reflectance 
changes due to snowmelt were more significant than differences among the dif- 
ferent types of melt-season sea ice. 

Also the spectral-hemispherical (plane) albedo of each measured arctic sur- 
face was computed. Comparing measured nadir reflectance to albedo for sea ice 
and snow-covered tundra shows albedo underestimated 5-40%, with the largest 
bias at wavelengths beyond 1 pm. For snow-free tundra, nadir reflectance under- 
estimates plane albedo by about 30-50%. 



1. Introduction 


The interaction of solar radiation with various arctic surfaces is an important 
component of the arctic climate. A measure of this interaction is the surface bidi- 
rectional reflectance function (BRF). The full angular distribution of BRF for a 
given surface (at a particular solar incidence angle) is known as the bidirectional 
reflectance distribution function (BRDF). Knowledge of the BRDF for a surface 
gives a complete description of its reflectance. Remote sensing instruments typi- 
cally have narrow fields of view, and measure only the reflected radiance at a 
single angle (thus only a small portion of the BRDF). Even the next generation 
satellites [King, 1999] will measure reflectance at a limited number of view an- 
gles. Thus proper interpretation of the remote sensing data will require knowl- 
edge of the BRDF of all surface types. 

The Cloud Absorption Radiometer (CAR) [King et al, 1986] flying aboard the 
University of Washington C-131A aircraft has provided a convenient and effi- 
cient means of obtaining BRF data. CAR BRF data over the past decade have 
been collected for several surface types. BRF data collected during the Kuwait 
Oil Fire Smoke Experiment (KOFSE) over an oil fire generated smoke layer is de- 
scribed by King [1992]. Tsay et al. [1998] describe BRF data gathered in Brazil over 
cerrado, dense forest, and smoke over dense forest, during the Smoke Clouds 
and Radiation - Brazil experiment (SCAR-B). Additional BRF data for the Saudi 
Arabian desert, forested wet land (Great Dismal Swamp), and ocean water con- 
taining sunglint over the Atlantic Ocean and Persian Gulf have been discussed 
by Soulen et al. [1999]. In this work we present BRF data of four common arctic 
surfaces: snow-covered sea ice, melt-season sea ice, snow-covered tundra, and 
tundra shortly after snowmelt. The snow-covered sea ice and tundra data were 
collected as part of the arctic Lead Experiment (LEADEX), and the melt-season 
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sea ice and tundra BRF data were collected during the Arctic Radiation Meas- 
urements in Column Atmosphere-Surface System (ARMCAS) experiment. 

Several studies of sea ice reflective characteristics have been reported [Gren- 
fell and Maykut, 1977; Grenfell and Perovich, 1984; Perovich et al., 1986; Allison et al, 
1993; and DeAbreu et al, 1995]. These sea ice studies, using typically ground- 
based instrumentation in mostly diffuse light conditions (overcast sky), discuss 
sea ice characteristics and their spatial and temporal effects on both spectral and 
wavelength-integrated albedo. Highlighted in particular are the changes that oc- 
cur with the onset of the melt-season. Prior to the melt-season (usually about late 
May in the arctic), sea ice albedo is mostly a function of the optical properties of 
the "dry" snow cover that overlies the sea ice. With the onset of the melt-season, 
however, the snow cover quickly begins to melt. Sea ice albedo then becomes a 
function of the sea ice characteristics such as brine and bubble concentrations, ice 
temperature, and ice thickness, coupled with the local meteorology (freeze-thaw 
cycles). The result is increased spatial variability and reduced reflectivity at all 
wavelengths (especially above about 1 pm where nearly all the incident solar ra- 
diation is absorbed). 

Direct application of these sea ice studies to remote sensing data (due to the 
characteristics of the remote sensing data) have heretofore been limited. To 
make the limited view angle remote sensing data more directly applicable, 
knowledge of the BRF at all azimuthal and polar viewing angles is needed. BRF 
of snow has been addressed by several studies [Dirmhirn and Eaton, 1975; Kuhn, 
1985; Steffen, 1987; Dozier et al, 1988; Hall et al, 1993; and Warren et al, 1998 (in- 
cludes a comprehensive review of factors affecting the BRDF of snow)], but BRF 
of sea ice has been limited. Perovich [1991] found that at azimuth angles 0° and 
90° from solar incidence, reflectance of a refrozen melt pond and young sea ice 
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0.2 m thick was relatively constant. Reflectance at zenith angles out to about 50° 
was relatively constant, but increased sharply for larger polar angles. Perovich 
[1994] reports that for sea ice progressing from cold dry snow-covered (pre-melt) 
conditions through various stages of melting, BRF data in the wavelength range 
0.4 — 1.0 pm at a polar angle of 30° is nearly azimuthally isotropic, but has a nota- 
ble specular peak at 60°. Both of these studies, however, were limited in wave- 
length range and number of polar angles measured. 

BRF studies of arctic tundra have also been limited. Recently, however, 
Vierling et al. [1997] and Eck et al. [1997] have reported results of ground-based 
BRF measurements taken during ARMCAS. They used the portable apparatus 
for instrument rapid acquisition of bidirectional observations of the land and at- 
mosphere (PARABOLA) to measure the BRDF of three types of arctic tundra: 
wet sedge (northernmost tundra), non-woody tussock tundra, and woody tus- 
sock tundra. The PARABOLA measures the nearly complete BRDF of three 
wavelength regions, but with an instantaneous field-of-view of 15°. Eck et al. 
[1997] report total albedo computed from the PARABOLA measurements to 
compare within 7.5% of pyranometer based measurements. Vierling et al. [1997] 
discuss how to use the BRDF measurements to differentiate between the three 
types of arctic tundra. 

In this paper the BRF data for each arctic surface are also used to compute 
surface hemispherical albedo. Surface hemispherical albedo represents the 
"true" albedo for a given surface. However, nadir reflectance (often the pa- 
rameter measured from an aircraft or satellite) is sometimes used as a surrogate 
[Barnsley et al, 1994]. Perovich [1994] reports that for sea ice evolving through the 
different melting stages, the spectral shape of nadir reflectance for wavelengths 
from 0.4 to 1.0 pm is similar to the spectral albedo shape. The ratio of these two 



Arnold et al.: Surface- Atmosphere Anisotropy of Sea Ice and Tundra 


4 


parameters was high (near 1.0) for the dry snow but, as melting progressed, the 
ratio decreased, especially at longer wavelengths. Similar comparisons have 
been done from our measurements and described below. 

In the next section, we quantitatively define surface bidirectional reflectance 
(as used in this work) and how hemispherical albedo is computed. In section 3, 
the instrumentation and data strategy are summarized, followed by discussion in 
section 4 of each measured surface spectral anisotropy and hemispherical albedo 
computation. A summary of this work is given in section 5. 


2. Definitions 

Following Tsay et al. [1998], we characterize the surface anisotropy from 
spectral measurements following the definition of bidirectional reflection func- 
tion (BRF, RJ found in van de Hulst [1980] and given by 




( 1 ) 


where 0 and 0 O are the viewing and illumination zenith angles respectively; <]> is 
the azimuth angle between the viewing and illumination directions; Po = cos 0 O ; 
F x is the collimated irradiance; and I x is the corresponding reflected radiance. 
The collimated irradiance F x is computed by weighting the solar flux at the top of 
the atmosphere with the spectral response function of each band, taking into ac- 
count the eccentricity of the earth's orbit around the sun. (Note: this definition 
does not account for atmospheric effects, but for most CAR channels this should 
be relatively small - see section 3 for additional details.) 

Similarly for the spectral albedo, r ^ radiance received by the detector from 
all viewing directions /solid angles are integrated (summed up for all IFOVs). 
Using Eq. (1), the spectral albedo is given by 
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where ( i , j) are the indices of (0, <|>) for discretization (M, N terms) respectively, 
over the hemisphere. 
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3. Instrumentation and data collection 

During LEADEX and ARMCAS, measurements of surface-atmosphere bidi- 
rectional reflectance were acquired by using a multispectral scanning radiometer 
- the Cloud Absorption Radiometer (CAR), which was designed and built at 
NASA Goddard Space Flight Center [cf. King et al., 1986 for details]. The CAR 
was mounted in the nose of the University of Washington's C-131 A research air- 
craft, providing radiometric measurements at 13 discrete wavelengths between 
0.3 and 2.3 pm. The first 7 channels are continuously recorded, whereas channels 
8 through 13 are located in a filter-wheel that can be locked at a particular chan- 
nel or manually or automatically rotated to measure a new wavelength interval 
after a preset number of scans. The wavelength location of each channel was se- 
lected to minimize atmospheric absorption effects (except for the UV-B channel 
where it is desired). The visible channel reflectance is weakly reduced by ozone 
absorption, and well as by aerosol and Rayleigh scattering. Aerosol effects 
should be small due to the relatively clean arctic air. The shortwave-infrared 
channel reflectance can be reduced by 2-8% due to water vapor and trace gas ab- 
sorption. 

Viewing with an LFOV of 1°, the CAR scans in a vertical plane on the right 
hand side of the aircraft from 5° before zenith to 5° past nadir (190° aperture), 
thereby permitting observations of both the zenith and nadir radiances to be ob- 
tained with as much as a 5° aircraft roll, an angle that is measured simultane- 
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ously with a gyroscope aboard the aircraft. Table 1 summarizes the CAR specifi- 
cations. 

The C-131A flew at an average speed of about 80 m s-1 in a circle with a 
minimum diameter of about 3 km, at a comfortable roll angle of 20° in about two 
minutes. At an altitude of 600 m above the targeted surface area and 1° IFOV, 
the pixel resolution is about 10 m at nadir and about 270 m at an 80° viewing an- 
gle from the CAR (see King [1992] and Tsay et al. [1998] for additional details). 

4. Results from observations 

The assumption of homogeneity of the targeted surface area is fundamental 
to properly reconstructing the bidirectional reflectance from CAR measurements. 
Thus targets were carefully selected to be as homogeneous as possible. During 
ARMCAS and LEADEX, eight spectral BRF data sets were obtained and ana- 
lyzed. The location of each BRF study on the north slope of Alaska and nearby 
Arctic Ocean is shown in figure 1, which further designates the precise location 
by assigning a letter from A-H. These eight studies have been categorized ac- 
cording to one of four types of surfaces measured: (i) snow-covered sea ice, (ii) 
melt-season sea ice, (iii) snow-covered tundra, or (iv) snow-free tundra. Two 
snow-covered sea ice BRF studies (A and B in figure 1) were acquired over the 
Beaufort Sea under nearly clear-sky conditions, from 20:50 to 20:57 and 21:15 to 
21:24 UTC on April 7, 1992. Three melt-season sea ice BRFs (C, D, and E) were 
also acquired over the Beaufort Sea under nearly clear-sky conditions, from 00:44 
to 00:55 UTC on June 6, 1995, and from 18:34 to 19:39 and 20:58 to 21:19 UTC on 
June 8, 1995. One snow-covered tundra BRF (F) was acquired just south of 
Prudhoe Bay, with nearly clear sky above (only some patchy thin cirrus), from 
23:38 to 23:43 UTC on April 7, 1992. Finally, two snow-free tundra cases (G and 
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H) were acquired south of Prudhoe Bay, with nearly clear sky above, from 06:37 
to 06:48 and from 06:57 to 07:09 UTC on June 12, 1995. 

4.1 Snow-covered sea ice (April 1995) 

Two three-orbit BRFs under clear sky conditions were conducted over snow- 
covered multi-year sea ice (locations A and B in figure 1). The solar zenith angle 
0 O for the two cases varied from 65.14°-65.26° and 65.22°-65.35°, respectively (cf. 
table 2). Figure 2a shows a representative photograph of snow-covered sea ice 
taken in the area where these BRF measurements were obtained. Due to the op- 
tically thick snow cover, the optical effects of the underlying sea ice is limited. 
Variations in surface reflectance are primarily due to shadowing caused by the 
relatively low sun angle and irregularities in the surface (pressure ridges). These 
effects are relatively small. The CAR filter-wheel channel was locked at 1.64 pm 
for 1 orbit in BRF A and for one orbit in B at 2.20 pm. 

Figure 3 contains the spectral BRF for each measured CAR channel over the 
snow-covered sea ice (since the BRF plots for A and B were similar, only those 
from A are shown here except for the 2.2 pm plot). In these and all other polar 
plots, zenith angle is represented as radial distance from the center, and azimuth 
angle as length of the arc on the respective zenith circle. The principal plane re- 
sides in the 0°-180° azimuthal plane with the sun located in the 0° azimuth direc- 
tion. With this definition, the upper half circle represents forward scattering and 
the lower half represents backscattering, with the anti-solar direction at 0 O « 65°, 
denoted by the white spot in the upper left plot of figure 3a. 

The BRF pattern in figure 3 shows reflectance minimized in the region 
around nadir and increasing going out from nadir, particularly in the solar and 
anti-solar directions. The strongest peak is in the forward scattering direction 
(specular reflection), while there is a diffuse and weaker peak in the backscatter- 
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ing direction. Interestingly, while the plot suggests a high degree of symmetry, 
reflectance mid-way between the solar and anti-solar directions is notably lower 
on the right than on the left in the dataset. The cause of this asymmetry is not 
clear. 

To investigate further the detailed angular distribution of BRF, figure 4a 
shows the spectral dependence of the reflection function along the principal 
plane. Except for the forward and backscattering regions near the limb, the plots 
for each band are otherwise fairly flat, representing a nearly Lambertian reflec- 
tance of snow-covered sea ice away from the specular reflectance direction. 

4.2 Melt-season sea ice (June 1995) 

The three melt-season sea ice cases were obtained in different geographical 
regions (locations C, D, and E in figure 1) to compare different types of sea ice. 
BRF C was 4 orbits, D was 35 orbits, and E was 10 orbits. Solar zenith angles 
varied from 54.51°-55.26°, 53.28°-57.69°, and 48.93M9.45 0 respectively, for C, D, 
and E. The filter-wheel was locked at 1.64 pm for 2, 13, and 4 orbits respectively, 
for C, D and E, and for 1, 9, and 4 orbits at 2.2 pm. A photo of the melt-season 
sea ice (from June 8) is shown in figure 2b. The BRF from C and E were both of 
shorefast ice but from different parts of the coastline. (Both are likely young sea 
ice and somewhat surprisingly had fewer open leads than D). According to 
flight notes, the surface in E showed no evidence of open leads or melt ponds at 
all (although a few ground spots were noted due to shallowness of the sea in that 
area). Much of the variability of the surface is likely due to melting and re- 
freezing and varying thickness of the ice. As compared to snow-covered sea ice, 
surface characteristics of melt-season sea ice show much more variability, par- 
ticularly at shorter wavelengths. 
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Data for BRFs C, D, and E are listed in table 2, and the BRF plot from D is 
presented in figure 5 (C and E are not shown but show similar features). The 
melt-season BRFs, like the snow-covered sea ice BRFs, show pronounced maxima 
in the forward scattering direction, and secondary maxima in the backscattering 
direction (though less apparent than for the snow-covered sea ice cases). Unlike 
the snow-covered sea ice, however, they show a distinct specular reflectance 
(sunglint) just to the right of the solar principal plane and at an angle from nadir 
equivalent to the solar zenith angle (slightly offset from the principal plane due 
to the aircraft pitch). The most important difference though is that the reflec- 
tance values of all the CAR bands for the melt-season cases are lower than for the 
snow-covered sea ice. This is most pronounced for the 1.6 and 2.1 pm channels. 
The melt-season sea ice surface for these channels is absorbing nearly all of the 
incident sunlight. The principal plane plots in figure 4b show reflectance of 1-2% 
in these bands, considerably lower than the snow-covered sea ice in figure 4a. 

4.3 Snow-covered tundra (April 1992) 

Following the two snow-covered sea ice BRFs (locations A and B), the C- 
131A flew a two-orbit BRF pattern over snow-covered tundra to the south of 
Prudhoe Bay (F in figure 1). The solar zenith angle varied from 65.14°-65.35° 
(table 2), and the filter- wheel was locked at 1.64 pm for one orbit. A photo of 
snow-covered tundra is shown in figure 2c, with the Trans-Alaska Pipeline 
clearly evident in this photo. Figure 6 illustrates the two-orbit average BRF for 
this case, which again shows maximum reflectance in the forward scattering di- 
rection with some evidence for a secondary backscattering maximum, though 
less apparent than for snow-covered sea ice. The cause of the asymmetry, the 
minima in reflectance on the left-hand side of the circle is not clear. Reflectance 
values are generally slightly less than the snow-covered sea ice cases on this day. 
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but are still quite high, even at 1.6 pm. A plot of the principal plane data is given 
in figure 7a. Snow covered tundra is quite Lambertian except for the forward 
scattering (specular reflectance) direction. 

4.4 Snow-free tundra (June 1995) 

Two BRF studies (5 orbits each) were conducted over the tundra south of 
Prudhoe Bay on June 12, 1995 (G and H in figure 1). These were "midnight sun" 
flights (solar zenith varied from 80.63°-81.23° and 82.00°-82.61°, respectively). 
The filter-wheel was locked at 1.64 pm for 3 and 2 orbits, respectively for G and 
H and for 1 and 2 orbits, respectively, for 2.1 pm. The tundra vegetation com- 
mon to this part of the north slope is senescent wet sedge (little photosynthesis 
since snow cover has just melted) [Vierling et ah, 1997]. Figure Id shows a close- 
up photograph of this type of tundra surface. Though not apparent in this pho- 
tograph the wet sedge was in many places partially submerged. The BRF plot in 
figure 8 (from G) depicts a fairly symmetric reflectance pattern for each case, 
with a fairly narrow maximum in the forward scattering direction (probably 
specular reflection due to the water in the tundra) and a broader (and higher for 
bands 0.87 pm and longer) backscattering pattern. The plot is similar in appear- 
ance to the Brazilian cerrado (savannah) case described by Tsay et al. [1998], 
though with a much lower sun angle here. Note the relatively low reflectance in 
the nadir direction at all wavelengths. The principal plane plot for G is shown in 
figure 7b. 

4.5 Spectral Albedo 

The spectral BRF measurements obtained by the CAR contain very fine an- 
gular resolution (i.e., M = 90 samples in 0 and N ~ 200 in <(>, cf. table 1) which pro- 
duces virtually no computational error in discretizing the integral. By using Eq. 



Arnold et al.: Surface- Atmosphere Anisotropy of Sea Ice and Tundra 


11 


(2), it is straightforward to compute the spectral-hemispherical albedo for the sea 
ice and tundra surfaces. The results for each of the sea ice and tundra cases are 
listed in tables 2 and 3, respectively. For each BRF case, columns 1 and 2 in each 
table lists the central wavelength and bandwidth characteristics of each CAR 
channel (note the change in channel 1 between 1992 and 1995). Column 3 lists 
the corresponding solar flux computed by weighting the solar flux at the top of 
the atmosphere with the spectral response function of each channel. The solar 
fluxes vary slightly from case to case due to the eccentricity of the earth's orbit 
about the sun. In column 4, the solar zenith angle varies less them 2° for all but 
one BRDF case (the 35 circle average of location D ~ 4.4°). Thus change in solar 
zenith angle does not contribute significantly to the variation in the computed 
spectral albedo (column 5). This variation is due largely to surface inhomogene- 
ity. Column 6 lists the average reflectance value and variation at nadir for each 
BRF study, for comparison to the spectral albedo. 

Generally for sea ice cases, the spectral albedo values in table 2 are consistent 
within a season, but show important differences when comparing the snow- 
covered sea ice to the melt-season sea ice, particularly for the longer wave- 
lengths. This can easily be seen in the plot of hemispherical albedo as a function 
of wavelength in figure 9a. The data suggest separate curves for the two sea ice 
types. The difference in the curves is about 10% at 0.67 pm, increasing substan- 
tially with wavelength becoming about a factor of 2 at 1.03 pm, to a factor of 6 or 
more at 1.64 pm. Also shown in table 2, the greater variability of the melt-season 
sea ice albedo at the shortest wavelength is indicative of a more variable surface 
(open leads, melt ponds, various types of bare ice, etc.) at this time of year. 

Nadir reflectance values as a function of wavelength are presented in figure 
9b for each of the sea ice cases, with the corresponding ratio of nadir reflectance 
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to albedo shown in figure 9c. Similar to the conclusion of Perovich, [1994], the 
plots of hemispherical albedo and nadir reflectance as a function of wavelength 
(figure 9a and 9b) have nearly the same spectral shape. Comparing the albedo of 
sea ice to the nadir reflectance (figure 9c), often used as a surrogate for surface 
albedo in climate models, yields two conclusions. First, at essentially all CAR 
wavelengths, the nadir snow-covered sea ice reflectance is a better surrogate for 
albedo than the melt-season sea ice albedo. Second, for both sea ice types, differ- 
ences between hemispherical albedo and nadir reflectance exceed 10% only for 
CAR wavelengths greater than about 1 pm (differences can be as much as 45%). 

Shown also in figures 9a and b is the effect on the data of applying an at- 
mospheric correction for Rayleigh scattering and ozone absorption (dashed 
lines). Perovich (1994) and many others have shown that the albedo and nadir 
values for ice and snow surfaces should continue to approach unity as the 
wavelength decreases. Clearly the uncorrected reflectance of CAR channels 1 
(0.47 or 0.50 pm) and 2 (0.67 pm) show just the opposite effect. Thus a simple 
correction for Rayleigh scattering and ozone absorption (aerosol absorption is 
assumed to be small in the clean arctic air) is applied to the data from BRFs A,D, 
F, and G and plotted as dashed lines in figures 9a and b. This results in a more 
realistic shape to the curves. Note, however, that this is only a first order atmos- 
pheric correction; a more thorough atmospheric correction is left for future work. 

Table 3 lists the spectral albedo and nadir reflectance values for the tundra 
cases. The snow-covered tundra albedos are slightly smaller than the snow- 
covered sea ice albedos at all wavelengths [cf. table 2]. Possible causes of this 
bias include age differences of the snow, temperature differences (degree of wet- 
ness of snow), contamination, and the effects due to underlying differences in 
surface reflectance (function of snow depth). Grenfell et al. [1984] discuss these 
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types of effects in further detail. Since no samples were collected of the snow in 
either case, it is not possible to determine the relative role of each of these effects. 
The difference in snow-covered tundra albedo as compared the snow-free tun- 
dra, however, is much more apparent. The snow-covered tundra albedo values 
are much higher than the snow-free tundra albedos at all visible wavelengths, 
but this difference decreases with increasing wavelength, and actually reverses at 
1.64 pm [cf. figure 9a]. 

Comparison of the hemispheric albedo of the snow-covered tundra to the 
measured nadir reflectance, similar to the sea ice cases, indicates only a small 
underestimate of the true albedo for wavelengths less than about 1 pm. Differ- 
ences can be as much as 40% for observations above 1 pm. For the snow-free 
tundra, however, the underestimate is much greater [cf. figure 9c]. The underes- 
timate ranges from about 30-50% and, unlike sea ice, the largest differences occur 
at the shortest wavelengths. 

5. Summary 

An important part of a successful effort to monitor the Earth's climate in the 
EOS era is the proper treatment of the sun's interaction with each type of Earth 
surface. To accurately interpret satellite measurements of top-of-the-atmosphere 
radiance, it is necessary to know the bidirectional reflectance function of each 
Earth surface. Thus the radiative effects of natural surfaces should be taken into 
account for remote sensing of aerosols and clouds, and also accurate surface 
hemispherical albedos can be determined for input to global climate models. 

During the LEADEX and ARMCAS field experiments, surface spectral bidi- 
rectional reflectance of sea ice and tundra for both snow-covered (pre-melt) and 
melt conditions were measured using the Cloud Absorption Radiometer (CAR). 
Snow-covered conditions for sea ice and tundra were monitored on April 7, 1992. 
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Sea ice in melt conditions was monitored on June 5 and 8, 1995, and snow free 
tundra on June 12, 1995. Despite contamination of the CAR data due to Rayleigh 
scattering and minor gaseous absorption (aerosol scattering was assumed to be 
small), results show distinct spectral characteristics for all types of surfaces pre- 
sented here, particularly the seasonal bidirectional reflectance changes in sea ice 
and tundra surfaces. In all cases sea ice reflectivity decreases with the onset of 
the melt-season, with the degree of difference increasing with wavelength. Re- 
flectivity of snow-covered tundra also decreased following snow-cover melt for 
CAR channels below 1.64 pm, but actually increased by a factor of two at 1.64 
pm. In addition, detailed measurements of the angular distribution of spectral 
reflectance were used to compute hemispherical albedos for each surface. For 
sea ice and snow-covered tundra, the CAR data indicate the surface spectral ani- 
sotropy only becomes significant at wavelengths greater than about 1 pm. Snow- 
free tundra, on the other hand, exhibits large surface spectral anisotropy at all 
wavelengths, with the greatest anisotropy at the shortest CAR wavelengths 
(~0.47 pm). Underestimates of the true albedo ranged from 30-50% when using 
nadir reflectance as a proxy for hemispherical albedo. 

The CAR measurements, combined with boundary layer aircraft measure- 
ments of aerosol physical properties and other co-located observations, form an 
unprecedented data set that is expected to greatly facilitate developing more re- 
alistic radiative transfer models in the future. All CAR data from both the 
ARMCAS and LEADEX field experiments have been processed to level-lb (cali- 
brated and geo-located radiance). Information on where and how to obtain CAR 
data from these and other field experiments (including all BRF polar plots), as 
well as browse images, a CAR data user's guide, and software for unpacking and 
interpreting CAR data, may be found on the World Wide Web at 
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http://ltpwww.gsfc.nasa.gov/CAR. 
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Table 1. Cloud Absorption Radiometer specifications during field campaigns. 


Platform 

University of Washington C-131 A aircraft 

Ground speed 

80 m s' 1 (nominal) 

Total field of view 

190° 

Instantaneous field of view 

17.5 mrad (1°) 

Pixels per scan line 

395 

Scan rate 

100 scan lines per minute 

Spectral channels 

13 (7 continuously sampled and 6 in filter wheel) 

Spectral range 

0.472 (0.503) - 2.303 pm 

Data system 

8 channels @10 bits 

Calibration 

Integrating sphere on the ground 
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Table 2. Computed spectral albedo and measured nadir reflectance for snow- 
covered sea ice (A-B) and melt-season sea ice (C-E). Uncertainties do 
not include possible errors in calibration. 


Central wave- 

Bandwidth 

Solar flux 





length (pm) 


(W m" 2 pm~l) Solar zenith angle 

Albedo 

Nadir reflectance 

(A) Snow-covered sea ice (April 7, 1992) 





0.503 

0.016 

1901.7 

65.14° 

- 65.26° 

0.737 ± 0.0040 

0.680 ± 0.0262 

0.673 

0.020 

1523.5 

65.14° 

- 65.26° 

0.770 ± 0.0034 

0.718 ± 0.0292 

0.754 

0.019 

1281.1 

65.14° 

- 65.26° 

0.818 ±0.0029 

0.776 ±0.0315 

0.866 

0.020 

947.85 

65.14° 

- 65.26° 

0.782 ±0.0012 

0.720 ±0.0286 

1.031 

0.020 

706.75 

65.14° 

- 65.26° 

0.710 ± 0.0036 

0.627 ±0.0264 

1.220 

0.021 

497.66 

65.14° 

- 65.26° 

0.506 ± 0.0060 

0.422 ± 0.0221 

1.270 

0.020 

410.52 

65.14° 

- 65.26° 

0.408 ± 0.0058 

0.343 ± 0.0187 

1.640 

0.041 

232.27 

wmm a 

- 65.26° 

0.196 

0.134 ± 0.0144 

(B) Snow-covered sea ice (April 7, 1992) 





0.503 

0.016 

1901.7 

65.22° 

- 65.35° 

0.739 ±0.0010 

0.686 ±0.0127 

0.673 

0.020 

1523.5 

65.22° 

- 65.35° 

0.774 ± 0.0016 

0.728 ± 0.0164 

0.754 

0.019 

1281.1 

65.22° 

- 65.35° 

0.823 ± 0.0020 

0.788 ±0.0176 

0.866 

0.020 

947.85 

65.22° 

- 65.35° 

0.784 ±0.0013 

0.730 ± 0.0168 

1.031 

0.020 

706.75 

65.22° 

- 65.35° 

0.704 ± 0.0029 

0.630 ± 0.0186 

1.220 

0.021 

497.66 

65.22° 

- 65.35° 

0.496 ± 0.0059 

0.421 ± 0.01% 

1.270 

0.020 

410.52 

65.22° 

- 65.35° 

0.399 ± 0.0056 

0.342 ± 0.0166 


0.040 

82.693 

65.22° 

- 65.27° 

0.199 

0.137 ± 0.0132 

(C) Melt-season sea ice (June 5, 1995) 





0.472 

0.021 

1962.9 

54.51° 

- 55.26° 

0.737 ± 0.0061 

0.696 ± 0.0813 

0.675 

0.020 

1471.7 

54.51° 

- 55.26° 

0.666 ± 0.0023 

0.626 ± 0.0852 

0.869 

0.022 

943.60 

54.51° 

- 55.26° 

0.581 ±0.0024 

0.519 ± 0.0772 

1.038 

0.020 

671.13 

54.51° 

- 55.26° 

0.379 ± 0.0023 

0.298 ± 0.0427 

1.219 

0.020 

475.07 

54.51° 

- 55.26° 

0.189 ± 0.0025 

0.130 ± 0.0194 

1.271 

0.021 

430.81 

54.51° 

- 55.26° 

0.175 ± 0.0021 

0.119 ± 0.0188 

1.643 

0.041 

229.64 

54.86° 

- 55.26° 

0.032 ±0.0022 

0.019 ± 0.0085 

2.207 

0.040 

72.250 



0.034 

0.020 ± 0.0132 

(D) Melt-season sea ice (June 8, 1995) 





0.472 

0.021 

1961.8 

53.28° 

- 57.69° 

0.734 ±0.0191 

0.689 ± 0.1367 

0.675 

0.020 

1471.0 

53.28° 

- 57.69° 

0.640 ±0.0192 

0.592 ± 0.1293 

0.869 

0.022 

943.11 

53.28° 

- 57.69° 

0.538 ±0.0178 

0.470 ±0.1157 

1.038 

0.020 

670.79 

53.28° 

- 57.69° 

0.341 ± 0.0128 

0.265 ± 0.0684 

1.219 

0.021 

474.83 

53.28° 

- 57.69° 

0.170 ± 0.0097 

0.117 ± 0.0305 

1.271 

0.020 

430.58 

53.28° 

- 57.69° 

0.156 ± 0.0092 

0.105 ± 0.0274 

1.643 

0.041 

229.52 

55.19° 

- 57.69° 

0.031 ± 0.0020 

0.016 ± 0.0089 

2.207 

0.040 

72.210 

54.41° 

- 56.72° 

0.031 ±0.0024 

0.019 ± 0.0128 

(E) Melt-season 

sea ice (June 8, 1995) 





0.472 

0.021 

1961.8 

48.93° 

- 49.45° 

0.707 ± 0.0080 

0.679 ± 0. 0996 

0.675 

0.020 

1471.0 

48.93° 

- 49.45° 

0.619 ± 0.0110 

0.592 ± 0.1100 

0.869 

0.022 

943.11 

48.93° 

- 49.45° 

0.511 ± 0.0105 

0.465 ± 0.1117 

1.038 

0.020 

670.79 

48.93° 

- 49.45° 

0.319 ± 0.0068 

0.259 ± 0.0615 

1.219 

0.021 

474.83 

48.93° 

- 49.45° 

0.154 ± 0.0039 

0.109 ±0.0251 

1.271 

0.020 

430.58 

48.93° 

- 49.45° 

0.145 ± 0.0040 

0.101 ± 0.0232 

1.643 

0.041 

229.52 

49.22° 

- 49.45° 

0.028 ± 0.0025 

0.016 ± 0.0075 

2.207 

0.040 

72.210 

48.97° 

- 49.12° 

0.025 ± 0.0013 

0.017 ± 0.0111 
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Table 3. Computed spectral albedo and measured nadir reflectance for snow- 
covered tundra (F) and snow-free tundra (G, H). Uncertainties do not 
include possible errors in calibration. 


Central wave- 
length (jam) 

Bandwidth 

(pm) 

Solar flux 
(W m" 2 jam' 

1) Solar zenith angle 

Albedo 

Nadir reflectance 

(F) Snow-covered tundra (April 7, 1992) 




0.503 

0.016 

1901.7 

65.14° - 65.35° 

0. 698 ±0.0141 

0.663 ± 0.0253 

0.673 

0.020 

1523.5 

65.14° - 65.35° 

0.725 ± 0.0189 

0.698 ± 0.0338 

0.754 

0.019 

1281.1 

65.14° - 65.35° 

0.770 ± 0.0214 

0.755 ± 0.0376 

0.866 

0.020 

947.85 

65.14° - 65.35° 

0.732 ± 0.0218 

0.699 ± 0.0395 

1.031 

0.020 

706.75 

65.14° - 65.35° 

0.659 ±0.0221 

0.606 ± 0.0403 

1.220 

0.021 

497.66 

65.14° - 65.35° 

0.462 ± 0.0181 

0.402 ± 0.0295 

1.270 

0.020 

410.52 

65.14° - 65.35° 

0.379 ±0.0146 

0.330 ± 0.0256 

1.640 

0.040 

232.27 

65.14° - 65.24° 

0.156 

0.104 ± 0.0126 

(G) Tundra - no snow (June 12, 1995) 




0.472 

0.021 

1962.9 

80.63° - 81.23° 

0.053 ± 0.0006 

0.027 ± 0.0070 

0.675 

0.020 

1471.7 

80.63° - 81.23° 

0.105 ±0.0017 

0.060 ±0.0117 

0.869 

0.020 

943.60 

80.63° - 81.23° 

0.266 ± 0.0041 

0.174 ± 0.0249 

1.038 

0.020 

671.13 

80.63° - 81.23° 

0.345 ±0.0054 

0.228 ± 0.0343 

1.219 

0.021 

475.07 

80.63° - 81.23° 

0.318 ± 0.0047 

0.214 ± 0.0335 

1.271 

0.020 

430.81 

80.63° - 81.23° 

0.271 ± 0.0038 

0.187 ± 0.0301 

1.643 

0.041 

229.64 

80.63° - 80.98° 

0.334 ± 0.0187 

0.209 ± 0.0428 

2.100 


92.280 

81.12° - 81.23° 

0.197 

0.121 ± 0.0370 

(H) Tundra - no snow (June 12, 1995) 




0.472 

0.021 

1962.9 

82.00° - 82.61° 

0.053 ±0.0026 

0.027 ± 0.0231 

0.675 

0.020 

1471.7 

82.00° - 82.61° 

0.096 ±0.0016 

0.056 ± 0.0292 

0.869 

0.020 

943.60 

82.00° - 82.61° 

0.260 ± 0.0007 

0.179 ± 0.0378 

1.038 

0.020 

671.13 

82.00° - 82.61° 

0.331 ± 0.0015 

0.231 ± 0.0479 

1.219 

0.021 

475.07 

82.00° - 82.61° 

0.298 ± 0.0025 

0.211 ± 0.0420 

1.271 

0.020 

430.81 

82.00° - 82.61° 

0.244 ± 0.0008 

0.176 ± 0.0367 

1.643 

0.041 

229.64 

82.00° - 82.26° 

0.328 ± 0.0054 

0.241 ± 0.0412 



92.280 


0.191 ± 0.0002 

0.128 ± 0.0461 
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FIGURE LEGENDS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6. 


Approximate location of each set of sea ice and tundra BRF measure- 
ments. 

Sample photos of each of the four surface types measured: (a) snow- 
covered sea ice; (b) melt-season sea ice; (c) snow-covered tundra 
(Trans- Alaska Pipeline and road in upper right corner of photo); and 
(d) tundra shortly after snow melt. Photos a-c were taken from the 
University of Washington C-131A aircraft, whereas photo d was taken 
from the ground. Photos by Michael D. King. 

Spectral measurements of surface-atmosphere bidirectional reflectance 
over snow-covered sea ice on April 7, 1992 during LEADEX. All bands 
are from BRF A (see figure 1) except i (2.2 pm) which is from BRF B. 
The location of the anti-solar point at 0 = 65° and <(> = 180° is indicated 
in figure 3a. 

Representative angular distribution of BRF in the principal plane for 
(a) the snow-covered and (b) melt-season sea ice surfaces. The capital 
letter after each channel label denotes from which BRF the principal 
plane plot was derived. 

Spectral measurements of surface-atmosphere bidirectional reflectance 
over melt-season sea ice. This case is BRF D from June 8, 1995 during 
ARMCAS, but is representative of the other two melt-season sea ice 
cases. The location of the anti-solar point at 0 = 55° and <j) = 180° is in- 
dicated in figure 5a. 

Spectral measurements of surface-atmosphere bidirectional reflectance 
over snow-covered tundra. This case is from BRF F from April 7, 1992 
during LEADEX. The location of the anti-solar point at 0 = 65° and <|) = 
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180° is indicated in figure 6a. 

Figure 7. Representative angular distribution of BRF in the principal plane for 
the (a) snow-covered and (b) snow-free tundra surfaces. The capital 
letter after each channel label denotes from which BRF the principal 
plane plot was derived. 

Figure 8. Spectral measurements of surface-atmosphere bidirectional reflectance 
over snow- free tundra. This case is BRF G on June 12, 1995 during 
ARMCAS. The location of the anti-solar point at 0 = 81° and <)> = 180° is 
indicated in figure 8a. 

Figure 9. Hemispherical albedo, nadir reflectance, and ratio of nadir reflectance 
to albedo plotted as a function of wavelength for BRFs A-H. The 
dashed lines represent the effect of adding a Rayleigh scattering and 
ozone absorption correction to the data (for clarity shown only for 
BRFs A, D, F, and G). 
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Figure 1. Approximate location of each set of sea ice and tundra BRF measure- 
ments. 
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Figure 2. Sample photos of each of the four surface types measured: (a) snow- 
covered sea ice; (b) melt-season sea ice; (c) snow-covered tundra 
(Trans- Alaska Pipeline and road in upper right corner of photo); and 
(d) tundra shortly after snow melt. Photos a-c were taken from the 
University of Washington C-131A aircraft, whereas photo d was taken 
from the ground. Photos by Michael D. King. 
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a) Snow-Covered Sea lce(0 o = 65") 




0.472 nm (D) 
0.675 nm (D) 
0.869 nm (D) 
1.038 pm (D) 
1.219 jim (D) 
1.271 Kim (D) 
1.643 fim (D) 
2.207 Kim (D) 


Figure 4. Representative angular distribution of BRF in the principal plane for 
(a) the snow-covered and (b) melt-season sea ice surfaces. The capital 
letter after each channel label denotes from which BRF the principal 
plane plot was derived. 
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Figure 7. Representative angular distribution of BRF in the principal plane for the 
(a) snow-covered and (b) snow-free tundra surfaces. The capital letter 
after each channel label denotes from which BRF the principal plane 
plot was derived. 
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Figure 9. Hemispherical albedo, nadir reflectance, and ratio of nadir reflectance 
to albedo plotted as a function of wavelength for BRFs A-H. The 
dashed lines represent the effect of adding a Rayleigh scattering and 
ozone absorption correction to the data (for clarity shown only for BRFs 
A, D, F, and G). 



